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Abstract 

This article studies a discrete geometric structure on triangulated manifolds and an 
associated curvature flow (combinatorial Yamabe flow). The associated evolution 
of curvature appears to be like a heat equation on graphs, but it can be shown 
to not satisfy the maximum principle. The notion of a parabolic-like operator is 
introduced as an operator which satisfies the maximum principle, but may not 
be parabolic in the usual sense of operators on graphs. A maximum principle is 
derived for the curvature of combinatorial Yamabe flow under certain assumptions 
on the triangulation, and hence the heat operator is shown to be parabolic-like. The 
maximum principle then allows a characterization of the curvature as well was a 
proof of long term existence of the flow. 
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1 Introduction 


In [9] we introduced the combinatorial Yamabe flow on three-dimensional 
piecewise linear complexes as an analogue of the smooth Yamabe flow (see 
[11], [20] for the Yamabe flow and [13] for a look at the Yamabe problem). 
The complexes are geometric in the sense that each Euclidean tetrahedron is 
given a metric structure by having the edge lengths determined by weights 
Tj dehned at each vertex i. The length of an edge {i,j} is dehned to be 
Tj -|- Tj] all such structures are called a conformal class as they are a conformal 
deformation of the triangulation where all edges are the same length. For 
each (nondegenerate) tetrahedron, we can think of the structure as coming 
from four mutually tangent spheres such that the centers are connected by 
the edges of the tetrahedron and the Rh sphere has radius rj. We showed that 
under sufficient long term existence conditions, the flow converges to constant 
curvature. 
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The three-dimensional combinatorial Yamabe flow was inspired by the work 
of Chow and Luo [3] (see also [14]). They looked at a combinatorial Ricci 
flow on two dimensional simplicial complexes and showed that the equation 
satisfies a maximum principle. The maximum principle, which says that the 
maximum of a solution decreases and the minimum increases, is one of the 
most useful concepts in the study of the heat equation and other parabolic 
partial differential equations. Maximum principle techniques have been used to 
great benefit in the smooth category. We are especially inspired by Hamilton’s 
work on the Ricci flow (see [12] and [2]). It is in general very difficult to prove 
a maximum principle, or even to prove short term existence of solutions, when 
you do not have a strictly parabolic equation. 

In this paper we investigate an analytic result about the flow which leads to 
a long term existence result for some structures. We And that the evolution 
of curvature admits a maximum principle under certain assumptions on the 
triangulation. It is especially interesting that we are able to derive a maxi¬ 
mum principle even in a situation where the evolution is not parabolic in the 
usual sense of graph Laplacians. We thus introduce the notion of parabolic- 
like operators which satisfy the maximum principle for a given function. We 
can then show that under sufficient assumptions the evolution of curvature is 
parabolic-like. 


2 Parabolic-like operators 


The weighted (unnormalized) Laplacian on a graph G = (V, E ), where V are 
the vertices and E are the edges, is defined as the operator 

(A/).= E “<U/i-/.) 

{tiles 


for each i G R, where the coefficients satisfy a^- = aji and a^- > 0 (see, for 
instance, [4]). The operator A takes functions on V to functions on V. The 
coefficients depend on the edge {i,j} (compare [5]). The symmetry condition 
is simply self-adjoint ness with respect to the Euclidean metric, so we can 
replace it with another self-adjointness condition. That is, we can define an 
inner product {f,g)i, = Y.ibifigi if we are given coefficients {hi} . An operator 
S is self-adjoint with respect to b if 

{Sf,g), = {f,Sg),. 

It is clear that symmetry corresponds to being self-adjoint with respect to the 
inner product determined by 6* = 1 for all i. We shall call {hi} a (positive 
definite) metric if hi >0 for all i. In order to match with the notation in [9], 
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we shall let So denote the vertices and iSi denote the edges. In later sections 
we will usually consider the inner product coming from the metric • 

We define a (discrete) parabolic operator on functions 

f:Sox[A,Z)^m, (1) 

where Sq is a discrete set of vertices and [A, Z) C M, as follows. First we shall 
call T the class of functions of the form (1). We write the evaluation of / at 
the point (i, t) as fi (t). 


Definition 1 An operator 
of the form 


P -T 


(Prkit) = ^- i: a,At)(f,{t)-f.(t)), 

where Oij : [A, Z) —M are self-adjoint with respect to some metrics {hi (t)}-^^^, 
is called parabolic if Oij (t) > 0 for all i,j G So and for all t G [A, Z). 

Parabolic operators are of the form ^ — A if A is dehned to be an appropriate 
Laplacian with weights. We note that it is easy to prove a maximum principle 
for these operators as follows: 


Proposition 2 If P is parabolic and f is a solution to Pf = 0, then f satisfies 

dfhi 


dt 

dfm 

dt 


(i)<0 

W>0 


where M,m E Sq such that 


fM (t) = max{/i (t)} 
fm{t) =min{/i(t)}. 

iGoo 


PROOF. Since 


and since umj > 


Pf = 0 we have, for a given f, 

^ = 5] aMj (t) {fj (t) - fM {t)) 

{M,j}esr 

0 and fj (t) < fM (t) for all i E So — {M} we see that 

%(i)<0. 
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The argument for fm {t) is similar. □ 


We may find operators that are of the form 

H ap (t) U'j {t) - fi {t)) 

{i,j}eSi 

but some coefficients are negative. The argument above does not work, but it 
is possible that a maximum principle still holds if the sum is positive when 
fi is minimal and the sum is negative when fi is maximal even though each 
term is not positive or negative respectively. This motivates our dehnition of 
parabolic-like operators as operators which satisfy the maximum principle for 
some function. 


Definition 3 An operator 

P 

of the form 

(•p/).(«) = ^- i: (AW-/.(*)), 


where aij : [A, Z) —M are self-adjoint with respect to some metrics {hi , 

is called parabolic-like for a function g if Pg = 0 implies 


dgM 

dt 

dg-m 

dt 


< 0 
> 0 


for all t G [A, Z). 


Parabolic-like operators formally look the same as parabolic operators, but 
some of the coefficients may be negative. We shall show that our discrete 
curvature flow equation is parabolic-like for the curvature function in a large 
subset of the domain. The hope is that this will be enough to prove pinching 
and convergence theorems. 


We note that the maximum and minimum may be done separately, dehning 
upper parabolie-like and lower parabolie-like in the obvious ways. In some 
situations we may be interested in only one of these. 


3 Combinatorial Yamabe flow 


Here we reintroduce the concepts of combinatorial Yamabe flow, based on 
the work on the combinatorial Ricci flow in two dimensions by Chow-Luo [3] 
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and the combinatorial scalar curvature by Cooper-Rivin [ 6 ]. Further details 
can be found in [9]. Recall that if 5 = {iSq, iSi, 1 S 2 , 1 S 3 } is a simplicial complex 
of dimension 3, where Si is the i-dimensional skeleton, we dehne the metric 
structure as a map 

r ; So^ ( 0 , 00 ) 

such that for every edge {i,j} G between vertices i and j, the length of the 
edge is iij = Ti+Vj. The set of all such metrics is called the conformal class since 
rescaling the r* will deform the structure to the metric structure with all edges 
the same length. We shall use T = {Tq, Ti, T 2 , T^} to denote the triangulation 
of one tetrahedron. Recall that in order for 4 positive numbers ri,rj,rk,re to 
dehne a nondegenerate tetrahedron, they must satisfy the Descartes inequality 


Qijkl 



Ti Tj Tk re 



- 2 



> 0 . 


We refer to Qijke as the nondegeneracy quadratic. Qijke is related both to the 
volume of the tetrahedron and the radius of the circumscripted sphere, i.e. 
the sphere which is tangent to all six edges of the tetrahedron. For a given 
tetrahedron, the existence of a circumscripted sphere is equivalent to being 
able to dehne the lengths by assigning weights r* to the vertices. 


The curvature Ki associated to a vertex i is dehned as 

~ 4:71 ^ ^ Oiijki 

{i,j,k,i}£S3 


where aijke is the solid angle vertex i in the tetrahedron {i, j, k,i} . The solid 
angle is the area of the triangle on the unit sphere cut out by the planes 
determined by {i, j, k} , , {i, k,i} where i is the center of the sphere. 

Note that aijke is symmetric in all permutations of the last three indices; when 
it is clear which tetrahedron we are working with, we will use the simplihed 
notation a*. The combinatorial Yamabe flow is dehned to be 


dn 

dt 


-KiVi 


Careful calculation shows that curvature satishes the following evolution 


dK 

-^= E (E - K^) + n^kji {Kk - Ki) + klujk {K, - Ki)]. 


This form is gotten by using the Schlahi formula (see, for instance, [17]) which 
can be written in the following way in this case: 


9oiijk£ doiijki dcXijki doiijki 
Ti—T^ -h -h Tk^:^ -h Te ^ = 0. 


dvi 


drj 


dvk 


drp 


( 2 ) 
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We computed the partial derivatives of the angles to be 


9cxij]f£ 

dvi 


Sr'jrlrj 


^PijkPijC-PiklVijkt 


2 111 

- 1 - 1 - 1 - 

Ti r-j Tk Ti 



1\ 

1 ^k , 

1 

1 

1 

H- 

1 + — 

- + - 

H- 

G VG rk 

reJ 

n 

VA 

A 

r£ 


H-(-1-1-j + (2rj + Vj + Tfc + ri) Qiju 

\ ^ j ^k / 


and 


( 3 ) 


9cXijk£ 

dvj 


.p‘T'y (- 

^PijkPijeVijke \ri \rj J r-j \ri r^J 



( 4 ) 


where Pijk is the perimeter of the triangle {i,j, k} and Vijke is the volume of 
the tetrahedron {i,j, k,£} . The coefficients flijke are 


^ijk£ 

and are thus easily computed to be 

1 


OcXijki 

\liikf = Ti 


dvi 


k^ijki 


4:rirjrlrj 
3 Pij k Pij t P^ij k£ 
1 

^r-k 


Geometrically, we hnd that 



( 5 ) 


^ kf. 




{i,j,k,e}eS3 


Tit 




where £*, is the area of the dual faces and the sum is over all k and where 
duality comes from assigning the geometric dual to a tetrahedron to be the 
center of the circumscripted sphere. The evolution of curvature can be written 
compactly as 


dKj 

dt 


AK, 


if we define the operator A as 


A/. 


4 E 



fi)- 


( 6 ) 
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Note that this operator looks like 


H Op- U'j - fi), 

which is a Laplacian on the graph (5o,5i) with weights aij, except that it is 
possible for to be negative. It can be argued that A is a discrete analogue 
of the Laplace-Beltrami operator. We shall see that since ^ijke are not always 
positive, the maximum principle is not ensured. In the next section we explore 
the maximum principle for this Laplacian. 


4 Parabolicity of ^ — A 


The operator ^ — A, where the Laplace-Beltrami operator A is dehned by 
(6), may not form a parabolic operator as dehned in Dehnition 1. This is 
because the coefficients ilijke = ^ay, in fact be negative. We can see this 
by using our explicit calculation for the coefficients in (5). If we choose, for 
instance, ri = r 2 = = 1 and = 1/5 then we see that the tetrahedron is 

not degenerate since 

Qi234 = 8 > 0 

and that one coefficient is negative. 


75Pi23T’i24W234 

This indicates that the maximum principle will not hold in general. However, 
there are several indications that this operator is a good operator and that it 
might be parabolic-like for the curvature in the sense of Dehnition 3. First is 
that the matrix is negative semidehnite, where the nullspace is spanned by 
the vector (ri, r 2 , rs, r 4 ). This is shown in Appendix A. It was originally stated 
in [6] but the proof there is incorrect. It is interesting to note that the incorrect 
proof tried to show that the coefficients of the Laplace-Beltrami operator are 
all positive, which we have shown to be false by a difficult calculation. Clearly 
the case when the coefficients are all positive is a large set. 


We pursue a diherent understanding of why the maximum principle should 
hold. The first evidence is that the Schlafli formula tells us that 


9o'j^ij£ d(y.£ijj^ 
- \-rj^ -hrfc ^ H-rr 


dvj 


dr-j 


dvi 


dvi 


0 


while the formula (3) for calculated in [9] shows that 


dcjjjke ^ Q 
dvi 


7 



for any tetrahedron. Together with the following lemma, we shall see that 
there are large restrictions on when coefficients can be negative. 

Lemma 4 IfQabcd < 0 then r^, n > min {r^, r^} . 


PROOF. Suppose Vtabcd < 0, then 


iri + i 

Ta \rb Tc 


+ 


rdJ n \ra 


+ 





< 0 . 


eliminating the denominators and regrouping terms we get 

TdTc [{rd + n) Tc + in + rj Td] + Ta (r^fc {Vd + Tc) - H {Vd - < 0 

so in particular we need 


TdTc {rd + rc) - n {rd -Vcf <0. 


Solving for we get 


rdrc {rd + rj ^ rdr^ 
{rd-rcf ~ Vd-r-c 


> min{rrf,rj 


since rd + rc> \rd — Tc| and max {rc, r^} > \rd — rd . Since the initial expres¬ 
sion is symmetric in a and b, we get r^ > min (r^^, rd too. □ 


Corollary 5 If {i, j, k, i} G S3 and r* = min (r*, r^, r*,, r^} then 

ktijke S: 0 and ^ 0 


or, equivalently, 


drj 


< 0 and 


dr-i 


> 0 . 


5 Monotonicity 


We need some way to relate the coefficients Ilijke and the curvatures. We 
attempt this by trying to prove that the K^s are monotonic as functions of the 
r’s. This will turn out to be true for the two cases of the double tetrahedron and 
the boundary of a 4-simplex, though not in general. We consider a tetrahedron 
determined by (ri, r 2 , rs, r 4 }. First we prove a lemma about the kinds of 
degeneracies that can develop. 




Proposition 6 If Qijke —^ 0 without any of the r* going to 0, then one solid 
angle goes to 2n and the others go to 0. The solid angle Oj which goes to 2n 
corresponds to ri being the minimum. 


PROOF. Rewrite QijM as 
Qijke 


1 

1 

1 


/I 1 1 1\ 

- 

H - 

H - 

H-) “2 


\ri 

rj 

rk 

rej ’ 

{rf rj ri rf J 


1 

1 

1 


1 1 1 

1 


- + - 

H - 

—) + 

— I-1- 

H- 

- — 


rk 

re 

ri J 

r-j \ri rk 

re 

rjj 


(7) 


1/1 1 1 1 \ 1 

— I-1-1-j -|- 

rk \ri rj re ru) re \ri 


+ — + 


Tk 


1 


If r* is the minimum, then 


Ti 

1 

Ti 

1 


rk 

1 


1 1 

-> 0 

re rj 

- - — > 0 

re rk 

— - - > 0 . 

rk re 


Hence if Qijke = 0 then 


1 

rk 


1 

re 


< 0 . 


ri 


Now we look at the partial derivative 


d 


J^Qijke — - 2 - 1 - 1 -— 


dr. 


rk re r* 


So if Qijke = 0, we can always increase to make the tetrahedron nondegen¬ 
erate. 

Now we can categorize the degenerations. Notice that by the formula for vol¬ 
ume, 

^-^ijk^^ije sin Pijke 


^ijke 


U. 




we must have that if the tetrahedron degenerates, sin fdijke —^ 0 for all dihedral 
angles. Hence jdijke goes to 0 or to tt. Since 


dpijke 

2rirjrir‘e 

1 

1 

-2^ 

{ 1 

dri 

^-Pijk-Pije ^ijke 

~T2 ‘ 
L ' k 

“ l2 

'e 

ri 

\ri rk 


( 1 1^ 

\ (2 

1 

l\ 



_l_- 

- + - 

H- ) 



VI 

1 \ri 

rk 

re) _ 



1 1 
+-+- 


ri re rk re 


2 + 
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(see [9]), if is the minimum, then < 0. When the tetrahedron hrst 

becomes degenerate, we can increase r* to become nondegenerate again. But 
this indicates that in this case jSijke would decrease, so if Pijke. = 0, Pijke. would 
become negative in a nondegenerate tetrahedron. This is a contradiction, so 
we cannot have jSijke = 0 in the limit. Hence 


Pijki Pikji Piijk ^ 


and aijki = 2n. Since aijke + cnjike + ctkije + oi^jk < 2n in any tetrahedron (see 
proof in [8]), we must have 


O^jiki Oikij£ (^tijk 0 - 


□ 


Now we can prove a monotonicity formula for angles in a given simplex as 
follows. 


Lemma 7 r* < Vj if and only if aijke > Oijike- 


PROOF. It is equivalent to prove that the strict inequality Vi < vj implies 
aijke > Oijike and that rj = r-j implies aijke = oijike- The second statement is 
clear. Since we are only looking at one tetrahedron, we can use a* instead of 
Oiijke without causing confusion. 

Consider the path u (s) = (ui (s), (T 2 (s), <73 (s), 0-4 (s)) dehned by < 7 ( 5 ) = 
(rj, Tj, (1 — s)rk + sri, svk + (1 — s) re) , where < r^. We can think of a as 
a function of four variables, where a {ri,rj,rk,r() = ctj. Let a^ {x,y, z,w) = 
a {z, X, y, w) and 04 {x, y, z,w) = a {w, x, y, z). 

Consider the function 

A= ( 8 ) 

isTo 

By the Schlafli formula, we hnd that 

dA = ^ ttj dvi 

isTo 


so 


dA 

ai 

UT i 

&^A daj dai 

dvidvj dvi dvj 
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More details can be found in [9]. Now, consider 

D (s) = {a (s)) = [as (a (s)) - 04 (a (s))] (r^ - r^). 
as 

The path is constructed so that D (1/2) = 0. Since the solid angles are between 
0 and 2n, D (s) < 2ti {re — r^) if the tetrahedron is nondegenerate for a (s). 
Consider the derivative 

-^D (s) = {a (s)) = [<9303 (a (s)) - 2dia^ {a (s)) + <94a4 (a (s))] (re - rkf • 

If the tetrahedron is nondegenerate, then -^D (s) < 0 since the Hessian of A 
is negative dehnite (see Appendix A) and none of the a* are equal to zero. 
Now, as we move along the path starting at s = 0, either the tetrahedron 
degenerates for some Sq G (0, 1/2) or the tetrahedron is nondegenerate up to 
s = 1/2. Suppose the first degeneracy is at s = sq. Then either 

as (so) = min {at (sq)} 

2 = 1,...,4 

or not. If it is the minimum, then at the degeneracy, 0:3 = 27i and the other 
angles are 0. This cannot happen, however, since then D (sq) = 27i (re — r ^), 
the maximum possible, but the derivative is negative, since -^D (sq) < 0 
implies that D (s) must have been larger than its maximum for some s < Sq, 
a contradiction. If as (s) is not the minimum, then at a degenerate point, 
^3 = «4 = 0. So there exists a first point Si G (0, 1/2] where = ae such 
that the tetrahedron is nondegenerate for s G [0, si). Hence D (si) = 0. Since 
the tetrahedron is nondegenerate on [0,si), we have -^D < 0 for s G [0,si). 
Together with H (si) = 0 we have D (s) > 0 for s G [0,Si). In particular, 
D(0) > 0, i.e. ttfc > ae- 

The reverse inequality must be true as well since the argument is symmet¬ 
ric. □ 


The lemma has the following interesting geometric consequence for a conformal 
tetrahedron which is an analogue of the fact that in any triangle the longest 
side is opposite the largest angle. The author does not know if this statement 
is true for a general tetrahedron. 

Corollary 8 For a conformal tetrahedron with metric structure {ri, r 2 , r 3 , r 4 }, 
the side with the largest area is opposite the largest solid angle, and the side 
with the smallest angle is opposite the smallest solid angle. 


PROOF. This follows from the fact that the angle atjke is opposite the side 
with area (r^ -|- -f re). □ 
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The lemma can be used to show monotonicity of curvature for small triangu¬ 
lations as follows. 


Corollary 9 In the case of the double tetrahedron, ri < rj if and only if 
K, < K,. 


PROOF. This follows from the fact that = 47r — 2aijke- □ 

Corollary 10 In the case of the boundary of a A-simplex, ri < rj if and only 
tfK, < K^. 


PROOF. Once again it is sufficient to show that r* < rj implies Ki < Kj 
and ri = rj implies Ki = Kj. The latter is trivial. Let’s number the vertices 
{ 1 ,..., 5 } . Suppose ri < r2. Then 


Ki — An — 01234 — ttl235 ~ «1245 ~ ^1345 
K2 = An — 02134 — CI2135 ~ CI2145 ~ C12345- 


We know by Lemma 7 that 02134 < 01234, 02135 < 01235, and 02145 < 01245. 
We thus need only show 02345 < 01345. Consider the path 


as = 


rir2 


(1 - s) r2 sri 


,0,0,r5 


We notice that that the nondegeneracy quadratic Q {a (s)) is a polynomial in 
s with highest term 



Since the quadratic Q {a (s)) is concave the minimum for s G [ 0 , 1 ] must occur 
at s = 0 or s = 1 . But since Q {a ( 0 )) = (^1345 > 0 and Q {a ( 1 )) = (^2345 > 0 , 
Q {a (s)) > 0 for all s G [ 0 , 1 ] . Furthermore 


fL 

ds 


a {a (s)) = ( 9 iO (a (s)) 


rir2 (r2 - ri) 

((1 - s)r2 + 


which is negative for all s G [ 0 , 1 ] by formula ( 3 ) since r2 > ri. Hence 
o (a (0)) > o (a (1)), or 01245 > a2i45- ^ 


Proving a similar statement about larger complexes would be much harder, 
since we cannot pair up angles which are in the same or bordering tetrahedra 
as we do here. We shall call this condition the monotonicity condition for a 
tetrahedron {i,j, k, £}: 


ri < rj if and only if Ki < Kj. 


(MC) 
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The condition is true for an open set of triangulations. Unfortunately, MC is 
not necessarily preserved under the flow. For instance, if r* = rj but Ki < K, 
then ^ = -Kifi > —KjVj = The mo notonicity will counteract some of 
the potential degenerations of the flow, and thus allows proof of the maximum 
principle and long term estimates. 


6 Proof of the maximum principle 


Suppose we have a complex such that each tetrahedron {i, j, k, 1} satishes the 
monotonicity condition MC. Assume that r* < Vj^Vk < r^. We shall hrst look 
at the minimum. By Corollary 5, klijki, klikje, kl^jk are all nonnegative, so since 
Ki is the minimum curvature among {i, j, /c,f'}, 

^ijk£ T ^ikj£ (,^k T ^iijk ^ b* 

Now for the entire triangulation, we get that if = miUjg^Q {Ki} , then 

^ {Kj T kl^nkji {Kk T kli()jk {^Ki Kjyl)\ 

{mj,fcy}G<S3 

must be a sum of nonnegative numbers since in any tetrahedron containing 
U = ^m, is the smallest weight. 

We now look at the maximum. We want to show that 


^iijk (^Ki K^ T kl^jik (^Kj K^ T kl^kij (^Kk K^^j ^ 0. 

This is certainly true if kleijk, klejik, klmj are all nonnegative since is the 
largest curvature. Again using Lemma 4 we see that if klabcd < 0 then b i, so 
^eijk > 0. We are then left with the case that both flejik, klikij are negative or 
only one is negative. First consider the case when both are negative. In this 
case it is sufficient to show that 

^£ijk T ^£jik T ^£kij ^ b 

since in this case we have 

^£ijk ^ k}i£jik kl^kij 

^iijk i,Ki Kf}) ^ (klfijik T kl(.kij) ruin {(^Kj Kf'j , Kf'j} 

— ^£jik (^Kj K(^ klikij i^Kk K(^ 


since 


{K, - Ke) < min {{K, - K ,), {Kk - K,)} 
< max {{K, - Ke ), {Kk - Ke)} 
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which is nonpositive and —flejik, —^ikij > 0. So 

^Hjk (-^2 ^tjik (-^j S~ ^ikij (-^fc ^ 0 

and we are done. The inequality ^lajk + ^ejik + ^mj > 0 follows from the 
Schlafli formula, since 

o I o I o 9^Hjk dcitijk dciUjk 

^^Ujk T ^^£jik T ^^ikij T j; ^'j T j; ^'k 


dvi 

0Oi£ij k 

dre 


dvj ^ ' dvk 


re- 


This is nonnegative by formula (3) for 

Now suppose that only one is negative, say ^ikij < 0, then similarly it is 
enough to show 

^iijk T ^£kij ^ 0 

since then 

^£ijk ^ ^ikij 

^£ijk ^£') — ^£kij i^k ^£') 

and flijik {Kj — Kg) < 0. We argue from our explicit calculations. The sum 
^£ijk T ^£kij Is equal to 


^rgr]rlrf /1 ^ 1 ^ 1V 1 ^ 1 ^ _ l' 

3P£ikP£ijVijk£ \re \rg rj Vk) 




1/1 1 1 

+ — + 




1/1 1 1 

— —^ - ^— 

Tk \re r-j r* 


Jj n 


which, when simplifying and multiplying by a positive factor, has the same 
sign as 


111 

{Pejk + ‘^Peij) - 1 " {P£jk + Ptij) - 1 " {Ptij + ‘^P£jk) - 

nrj TiTk rgTi 

+ {‘^Pljk + P£ij) - 1 " {P£ij + P£jk) - 1 " {‘^P£ij + P£jk) - 

TjTk rgTj rgVk 

1 . 1 „ 1 


P£ij 2 i.P£jk T Piij) 


_p _ 

,2 ^2 ■ 


This is greater than or equal to 


P£ijQijk£ + Pejk 


JPj rj 


“2 I + {P£jk — P£ij) ' - 


nVk r% 
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which is nonnegative since r* < rj and {i,j, k,i} is nondegenerate. 
Thus we have proven the following: 


Theorem 11 On a complex with a metric structure which satisfies the mono¬ 
tonicity condition MC, if Km is the minimum curvature, Km is the maximum 
curvature, and we satisfy the combinatorial Yamabe flow 

-—r-i = -Kiri 

dt 

then 

—Km > 0 
dt 

and 

^Km < 0 , 

dt 

i.e. the combinatorial Yamabe flow is parabolic-like for K. 

Corollary 12 On the double tetrahedron and the boundary of a f-simplex the 
combinatorial Yamabe flow is parabolic-like for the curvature function K, i.e. 

d 

< 0 
dt 
d 

— Km > 0 . 

dt 


PROOF. This follows from Corollaries 9 and 10 which say that the mono¬ 
tonicity condition is satisfied in these cases. □ 


The maximum principle has been used by Hamilton and others to prove many 
pinching results for geometric evolution equations (e.g. [10], [12]). The most 
basic use is to show preservation of positive or negative curvature, which is an 
easy corollary. 

Corollary 13 On a complex with a metric structure which satisfies the mono¬ 
tonicity condition MC fort G [0,T), then nonnegative curvature is preserved, 
i.e. if Ki > 0 for all i for t = 0, then Ki > 0 for all i for all t G [0,T). 
Similarly, nonpositive curvature is preserved. 


PROOF. If > 0 for all i, then in particular the minimum is nonnegative 
and increasing. □ 
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7 Long term existence 


The monotonicity condition (MC) also gives us a way to show long term 
existence. The maximum principle will allow us to bound the growth or decay 
of the lengths rj. 


In order to show long term existence, we must show that Qijki > 0 for every 
{i, j, k, i} G ^3 and that r* does not go to zero or inhnity in hnite time for any 
f G iSq. Since we are only working with one tetrahedron, we can use Q = Qijki 
without fear of confusion. We calculate 


dQ 2/1 1 1 

7 ^ —- 2 - 1 - 1 - 

rf \rj Tk Ti r-i^ 




(9) 


and hence using formula (7) we see 


on (dQ dQ dQ dQ ^ 

+ —Tj + -^Tk + -^re 
\ uTi uTj ur^ UT£ j 


If Q = 0, then we have 


dQ dQ dQ dQ 
—n + —r-j + —Tk + —ri = 0. 


dr. 


dvi 


dvk 


drp 


( 10 ) 


Write the evolution of Q as 


dQ dQ dQ dQ dQ 

~Tr = - —KkVk - -^KeVe 

dt OTI UTj OT ^ OTI 

(ir, - Kj - ^Tk (Kk - Kj - {Kp - Kj 


dr 


drp 


using (10). Notice that if r* is the minimum, then ^ < 0 for j ^ i hy (9). If 
the monotonicity condition MC is satished, then we must have 


dQ 

dt 


> 0 


if Q = 0 since Kj > Ki and hence the tetrahedra do not degenerate. 


Now we can show that the maximum principle gives bounds on growth and 
decay of the r*. 

Proposition 14 If < Q and > 0 and then there is a constant C 
such that 

ri (0) < ri {t) < ri (0) e^f 
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PROOF. Let C (t) = max {Km (t) , —Km (t)} > 0. Then 

-Cvi < -KiTi < Cvi 

for each t. Since < 0 and > 0, we must have C (t) < C (0). So if we 
look at the evolution 

uT ■ 

-C^(0)r,<^<C(0)r, 

we get 

n (0) < n {t) < n (o) 

□ 


Thus the solution exists for all time. Convergence to constant curvature now 
follows from [9]. 


8 Further remarks 


In this paper we have seen two large sets of possible metric structures within 
which the maximum principle holds: the set where the coefficients are 
positive and the set where the monotonicity condition MC is satished. Un¬ 
fortunately, neither of these conditions is obviously preserved by the flow. It 
would be highly desirable to hnd a set which is preserved by the flow within 
which the curvature satisfies the maximum principle. 

Numerical data suggests that the maximum principle holds in much greater 
generality, even for large simplicial complexes that do not satisfy monotonicity. 
Numerical simulation of the flow requires a true simplicial complex; a CW 
decomposition will not work because there are not enough vertices to allow the 
different tetrahedra in the complex to evolve independently. Thus the current 
numerical work has been limited to certain small triangulations (fewer than 
15 vertices) of the 3-sphere, the direct product of the 2-sphere with the circle, 
the twisted product of the 2-sphere with the circle, and the 3-torus. Some of 
the small triangulations are due to F.H. Lutz (see [16] and [15]). The condition 
of monotonicity is not particularly well understood for large complexes either, 
though it is known not to hold in general even for triangulations of S^. 

The maximum principle is closely connected to the fact that the operator 
A is negative semi-dehnite in the smooth case, but it is not clear that the 
dehnition of maximum principle for graph Laplacians which we use here is 
the right maximum principle to correspond to the dehniteness. This may also 
be related to the fact that the principle eigenfunction in the smooth case 
is positive, while in the discrete case here the principle eigenvector usually 
will not have all positive (or all negative) entries. Perhaps a kind of ‘refined 
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maximum principle’ is needed, as in [1] and [18]. Also integral type maximum 
principles have been successful in studying discrete Laplacians as in [7]. 
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Appendix A 


In this appendix we prove that the four by four symmetric matrix has 

three negative eigenvalues and one zero eigenvalue. The proof in [6] is incorrect; 
Igor Rivin has given a new proof in [19]. We do not understand part of Rivin’s 
proof, but we can complete Rivin’s arguments by a calculation given below. 
The complete argument is given here. 

The zero eigenvalue comes from the vector (ri, r 2 , rs, r^) by the Schlafli formula 
(2). Recall that the (n — l)-dimensional minor M {i,j) of a matrix M is the 
matrix with the fth row and the jth column removed. Let A be the matrix 
of partial derivatives, so Aij = We take as the domain of A the set of 

(ri, r 2 , rs, Vi) such that the associated tetrahedron is nondegenerate, i.e. r* > 0 
for f = 1,..., 4 and Qijke > 0. 

Proposition 15 The minors A {i,j) have determinant 
detA(i,j) = (-!)*«+■ p 

and hence is nonzero if the tetrahedron {i,j,k,i} is nondegenerate. 


PROOF. We can do a rather lengthy calculation using (3) and (4). Note that 
we need only compute the minors A{i,j) where i ^ j since detA = 0. The 
minors A{i,i) on the diagonal are slightly more difficult to calculate because 
there are 3 entries of the more complicated form ^ instead of only 2 for the 
off-diagonal A {i,j) where i ^ j. □ 

Corollary 16 The matrix A is negative semidefinite, rank 3, and the nullspace 
is the span of the vector (ri, r 2 , rs, r 4 ). 
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PROOF. The rank follows from Proposition 15. The nullspace condition fol¬ 
lows from the Schlafli formula (2). Since the domain is connected and the rank 
is always 3, the eigenvalues must always have the same sign. We need only 
compute the matrix at one point, say rj = 1 for i = 1 ,..., 4. The matrix A at 
this point is easily computed to be 

-3111 
1 1-311 

3^2 11-31 ’ 

111-3 

which has eigenvalues 0, — 
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